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Abstract — We consider the problem of Joint Routing, Schedul- 
ing and Power-control (JRSP) problem for multihop wireless 
networks (MHWN) with multiple antennas. We extend the 
problem and a (sub-optimal) heuristic solution method for JRSP 
in MHWN with single antennas. We present an iterative scheme 
to calculate link capacities(achievable rates) in the interference 
environment of the network using SINR model. We then present 
the algorithm for solving the JRSP problem. This completes a 
feasible system model for MHWN when nodes have multiple 
antennas. We show that the gain we achieve by using multiple 
antennas in the network is linear both in optimal performance 
as well as heuristic algorithmic performance. 

I. Introduction 

Multihop wireless networks are essential for ubiquitous 
computation and communication. Currently there are many 
experimental setups of multihop wireless networks around 
the world. Ad-hoc wireless networks and sensor networks 
are also examples of multihop wireless networks where it 
is necessary to employ multiple wireless hops for even the 
connectivity of the nodes deployed in a particular area. How- 
ever, multiple wireless hops pose many new challenges in a 
network design. But recent studies have shown that some of 
these challenges can be converted into opportunities by careful 
network design. Thus new communication paradigms, e.g., 
opportunistic scheduling, cooperative communication, network 
coding and multiple antennas have been developed in recent 
years. Exploiting these techniques together in a multihop setup 
to optimize the system performance is very challenging. In 
this paper we will concentrate on designing multihop networks 
with multiple antennas at each node. 

In multihop wireless networks, unlike in the wire-line 
networks, even the concept of a link between two nodes 
is not properly defined. Shadowing can, irrespective of the 
distance between two nodes, cause the channel to be very 
bad. Even when there is no shadowing, we may or may 
not be able to communicate directly between two nodes 
depending on the transmit power used and also if there are 
other users transmitting in the neighborhood. Thus, topology 
of the network, transmit power, link scheduling and routing are 
all interrelated and for optimal performance one may need to 



jointly optimize power, scheduling and routing. This problem 
is computationally intractable and not scalable even for a 
centraHzed algorithm |[l|, (2), Q. 

Employing multiple antennas at a transmitter and/or at 
a receiver can provide transmit diversity, receive diversity, 
increase the capacity of the link and reduce BER. Thus, in 
wireless networks where bandwidth is scarce, it is important 
to employ multiple antennas wherever feasible. This increases 
the degrees of freedom one can exploit to improve the system 
performance. Thus even for multihop network it is desirable 
to have multiple antennas at different nodes. However, as 
mentioned above, even with single antennas, jointly optimizing 
power, routing and scheduling in a multihop wireless network 
is very computationally intensive, with multiple antennas the 
problem becomes extremely challenging. At the moment these 
are very few studies available for this system. We address this 
problem in this report. 

II. Related Work 

In wireless networks, it is well known that the traditional 
layers of the communication network cannot be considered 
in isolation. Many authors have proposed joint approaches to 
issues such as power control, scheduling, and routing (T\, (T\, 
|4 |, |5]. But there has been little effort in joint optimization of 
power, link scheduling and routing. One of the early papers in 
this direction is |T|. In fl], the techniques of (V\ are applied 
to a different scenario where the nodes use Energy Harvesting 
Sensors. Both ||Tj and ||2) use Linear Programming in their 
system models by using the technique of Column Generation. 

A multihop wireless network with single antennas is 
studied and solved in a similar way claiming to achieve 
the solutions for much larger networks than what has been 
possible thus far in |j9J. 

III. System Model And The Problem Formulation 

We have a network of half duplex nodes spread out in space 
communicating with each other using the common wireless 
medium. Since the network uses a wireless channel of broad- 
casting nature, ideally every node can hear the transmissions 
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from every other node. But in practice, a node usually will 
have a coverage area hence being allowed to transmit to or 
receive from a limited number of other nodes. This effectively 
forms a Graph and is the starting point for the problem. But 
it can be noted that interference continues to exist between 
any pair of nodes, irrespective of the graph representation. 
Further, the nodes may be restricted in terms of transmission 
or reception towards other nodes(neighbors) in the network. 
This is required for a tractable model in terms of formulating 
and solving the problem. We term it as transmission model of 
the system. 

In the given network, a subset of nodes called sources 
would like to transmit to their destinations, which is another 
subset of nodes. It in general needs multiple wireless hops to 
reach the destination. Our objective is to provide a fair (to be 
defined later) data transmission rate to all those sources via 
multiple intermediate nodes. 

A. System Model 

We have a network being represented as a fully connected 
directed graph with no self-loops, denoted as Q{J\f,£). Here, 
Af = {1,2, . . . , N} is the set of nodes and £ = {1, 2, . . . , L} 
is the set of directed links in the network. The set of source 
and destination pairs and the transmissions among them will be 
called as^ows. The set of flows is denoted by = 1, 2, . . . , F. 
Each flow f ^ J- will be associated with a source node and a 
destination node denoted by src(f) and dest{f). It is allowed 
to have common source and/or destination nodes among flows. 
Each f E T will also have a set of rate demands denoted 
by {di,i e J-}, while the system may decide to allow rates 
equal to {ri,i £ J^}, depending on the availability of network 
resources. The rates allowed may be lesser or greater or equal 
to the demanded rates. We define fairness as the least fraction 
in 



1 2^ I and denote it by A. 



All the nodes use multiple but equal number of antennas a 
for transmission or reception. The total power consumption of 
a node for transmission can only be a value chosen from a set 
of K powerlevels (including zero), in any time slot. We restrict 
all the antennas to be used towards a single node at any time 
instant. That means, we use point to point transmission model 
or the beam forming model of transmission in the network. 
Hence, a node involves in at most one active link, at any point 
of time. 

Though this model is a simple extension of the model used 
in |[T] and I^J for single antennas, it is not obvious to formulate 
and solve the problem for multiple antennas. Further, one 
cannot simply state the MIMO gain in the final throughput, 
just by observing that each link gains in its capacity, because 
of multiple hops being involved in transmissions. 

Given the transmission model of the network, not every 
subset of C can be activated simultaneously. We define a mode, 
as a valid subset of links being activated simultaneously along 
with the amounts of powers used on the links. We represent it 
by an i X 1 vector of powers used on each link. A set of modes 



given a schedule, achieves the link scheduling and routing 
simultaneously. If such scheduling takes the power availability 
on an average at each node into account, it becomes a solution 
to the JRSP problem. 

The fading channel gains are assumed to be known during 
the analysis being a centralized setup and are given by set of 
a X a matrices Hij{\/i,j e JV), which represent the channel 
matrices between nodes i and j. We also use the notation 
H'i^{yi,m e £) to denote the channel between links l,m 
meaning, the channel between source node of link I and 



destination node of link m. Note that H',, 



Hij, where 



i = src{l) and j — dest{m). Each scalar in the matrices is 
taken as a circularly Gaussian complex random variable. In 
|1| and |2], these were calculated by using a simple path loss 
model for a set of nodes randomly spread out over an area, 
based on the euclidean distances between them. 

Given the power spent on a link (in any mode), we retain 
the distribution of powers over the set of antennas at nodes 
as a degree of freedom for nodes (and hence is a variable 
during link capacity calculations). Hence they use waterfilling, 
whenever required. The way we compute link capacities is 
presented after the mathematical statement of the problem. 

B. The Mathematical Formulation 



max A = ( min s — > ) (la) 
(fairness objective) 

subject to : 

AX = Y , (lb) 
(Flow Conservation) 
X -l < C -a, (Ic) 
(Avg. Link Capacities) 
p. a < P^'^a ^ (Id) 

(Avg. Power Control) 
a • 1 = 1 , (le) 
(Consistency of scheduling) 

Pn > 0, a,n > 0, Vn e A/", Vm e X (If) 

(non-negativity of variables) 

The problem setup used in 11] and |2] is quite general 
and the problem formulation under multiple antennas can be 
made similar. Except that the terms used in it come from a 
different network scenario in multiple antennas. Its a conven- 
tional Multicommodity Flow problem from network theory and 
operations research literature, which is modified and exploited 
to suit the JRSP problem. The problem's objective is to find 
the solution to the JRSP problem, that gives maximum /a/rneii 
among the users. The mathematical statement of JRSP problem 
contains the well known flow conservation constraints, link 
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capacity constraints and average power constraints. This forms 
a simple linear programming problem (LPP) as given in ([T]). 

In the problem statement ([T]i, 'A' refers to the node versus 
link representation of the graph which takes the form of an 
N X L matrix with elements Uij taking signed binary values 

0, ±1, as below: 

{0, if node i is not a part of link j 
+1, if node i is the source to link j, 
— 1, if node i is the dest to link j. 

X refers to an L x matrix with values representing 
the average effective rate of transmission chosen on link I 
corresponding to flow /. 

r refers to an x F matrix, which contains the net 
outwards data rate at each node corresponding to each flow. 

1. e. 

{0, if node i is not a source or dest. of flow j 
+rj, if node i is the source to flow j, 
—Tj, if node i is the destination to flow j. 

Ai — {1, 2, M} refers to the set of all possible modes 
in the network. Each mode m G M will take a representation 
of an L dimensional column vector. The modes set includes 
an idle state (i.e. a 0) as well. 

a represents the vector of time fractions given to each of 
the modes. 

C is an M X L vector, where an element C[" represents the 
capacity of a link I G C in mode m e A4. i.e., each column of 
C is a vector of same dimension as a mode, but link powers 
are replaced with link capacities. 

P refers to an TV x M matrix called power profile matrix 
with values P™ representing the amount of power spent by 
node n in mode to e A^. It corresponds to the power spent 
on the single active link in the neighborhood of a node as 
a transmitter, at any time slot. P""^ is an x 1 vector, 
representing the power availability at each node, on an average. 

Note that the variables in the problem are ri(Vz E T),X 
and a while rest all(A, P, C and p''"^) are constants once the 
network and set of all modes A4 are fixed. Also note that, 
its a pure linear program since the objective can alternately 
be written as a set of linear inequalities introducing one new 
variable. 

With the set of assumptions we made in the system model, 
given a network(A) with average power values(P''^^), this 
problem statement becomes universally applicable for any 
number of antennas(a) on nodes, given the set of modes A4 
and their capacities C. But with multiple antennas, we see 
following issues. 

1) The link capacity calculations are interdependent (a joint 
optimization problem) due to interference. We would have 
a capacity region of points being vectors of link capacity 
elements. 

2) The capacity region of such a vector interference channel 



is still an open problem. 

3) We need a computationally feasible method to decide 
link capacities in a mode (columns of C), atleast as a good 
achievable rate point. Because, we have to calculate the same 
for a large number of modes. 

C. Link Capacity computation 

Let TO = (toi, TO2, . . . , tol)-^ be the given vector of L 
powers spent on links. To compute the link capacities, one 
should solve the problem of waterfilling on each of the links. 
But, the effective noise that each user sees is a sum of AWG 
noise and the interference from other users. Assuming each 
source emits independent Gaussian symbols, we need the 
covariance matrix of the effective noise, for which we need 
covariance matrices of each transmitter in the network. This is 
clearly not known at the time of computation of link capacities. 
Because, link capacity computation itself is water filing i.e. 
an optimization decision over all transmit covariance matrices 
possible at all the links. 

We propose the following iterative scheme to greedily 
calculate the capacities of links in the network which are 
active simultaneously. We continue the iterations till we see 
the convergence in sum rate of all links in the network. The 
sum rate increases monotonically after each iteration. 

We use the following notation. Ci is the i*'* link capacity. 
Ki is the transmit covariance matrix used on link i while 
TOi is the power spent on link i. sumrate corresponds to 
the sum of all link capacities active in the network. The 
function "waterfill(.)" performs waterfilling on any specified 
link and outputs two values namely link capacity and optimum 
transmit covariance matrix to achieve that capacity under the 
interference value at that instant. It takes three arguments 
namely link index at which waterfilling is to be performed, set 
of all transmit covariance matrices and channels(i7|^ ) between 
link i and every other link j in the network, (i.e., source of 
a link to the destination of other link — H'^j — Hpq where 
p — src{i) and q — dest{j)). 



Algorithm 1 Iterative Waterfilling algorithm for finding link- 
capacities in a point-to-point mode 

initialize sumrate — — oo, C; =0, Ki — (nii/a) * / 
Vz = {l,2,...,L} 
while I J2i Ci — sumrate] > e do 
sumrate = C'i 
for i = 1 to L do 

if link i is active then 

[K,,a] = waterfiU(i; a^; {if,, i/^.Vj £ £}) 
end if 
end for 
end while 



D. Complexity Issues 

The number of constraints in the problem is {NF + L + 
F +1) while number of variables is {F + LF + M). If N is 
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taken as the size of the network, N < L < N{N — 1) since 
graph is assumed to be directed and fully connected, having 
no self loops. Further F < N{N - 1). 

Though one can derive an upper bound on M as (a 
super exponential number) by arguing that M is less than 
all possible choices of powerlevels among L links, it can be 
seen via bruteforce enumerating simulations that the number 
of possible modes in a network of N nodes is exponential in N 
(see fig[T]). In anycase, the number of variables in the problem 
grows enormously with N and is much larger compared to 
number of constraints in the network. Typically, this number 
turns out to be more than a milhon for a network as small as 
1 1 nodes and 4 powerlevels. 

P2P 
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Fig. 1. A typical number of modes(M) as a function of number of nodes(N) 

Hence, its a prohibitively large problem to be solved 
directly for networks more than roughly 10 nodes. Note 
that, we have ignored the calculations required to get link 
capacities. 

IV. The Heuristic Solution 

One definitely needs to have a solution procedure, that is 
feasible for solving the network, upto atleast reasonable size 
of the network. [1| proposes a heuristic column generation 
method to solve this problem and is claimed to solve networks 
upto approximately 30 nodes. Q applies the similar thing with 
slight variation, on the networks of energy harvesting sensor 
nodes. We see that we can still solve the problem upto almost 
the same level, even when the network has multiple antennas 
and hence the complexity is much larger. 

The solution method proceeds as follows. First it can be 
seen by using the analysis of number of basic variables that, 
no more than N + L + 1 number of modes are necessary 
to obtain the optimal basic feasible solution. Then as per the 
column generation procedure, solution starts by considering a 
smaller problem called Master Problem with just iV + i + 1 
randomly chosen mode variables(a') from a, rest all being the 
same to the original optimization problem in ([T]i. 

It also considers a Sub Problem that chooses a new mode 
called a good mode from the remaining modes, considering 



which we would improve the performance of the master 

problem. 

Sub Problem: 

max e{m) ^ u^C„, - v^P"' - (3 

meM\M' 

subject to : 
eim) > 0. 

Such a mode will replace an already chosen mode variable 
which is given zero scheduling time in the optimal solution 
of Master problem. This starts next iteration, in which master 
problem attains an improved solution. 

The same after iterating for sufficient number of iterations, 
converges to a (sub)optimal solution and an (sub)optimal 
choice of set of iV + L + 1 modes from Ai. 

In this procedure, sub-problem is the stage, which requires 
an optimum choice of variable by searching exhaustively 
among auxiliary function values evaluated at A/— (N + L + l) 
modes, which is again computationally infeasible. This can be 
avoided by using a heuristic algorithm for sub-problem, from 
IT] and fT\. This comes at the cost of sub-optimality to the 
solution. But it was shown via simulation examples that, the 
solution that we see via this heuristic, is close to the optimal 
solution. 

The heuristic algorithm for MHWN with MIMO extended 
from its single antenna version is presented below as Algo- 
rithm |2j in a ready to implement algorithmic form. In this, 
nextlevel{x) is the notation used for denoting the power value 
which is least among the power values that are greater than 
X from power level set. And inter ferers{l) is the notation 
used to denote the set of all the neighbouring links to link I 
that share a node with link I. 



Algorithm 2 The Heuristic Algorithm to obtain an efficient 
sub-optimal solution to the sub-problem 

Initialize mode vector goodmode = 0, 
lasttheta = 9(goodmode), AllowedSet = C 

while AllowedSet ^ do 
for I = 1 to L do 

if i e AllowedSet then 
m = goodmode 

nii = nextlevelirrii) 
di = e{m) 
end if 
end for 

Let besttheta = max^ di & 
Let the bestlink I = arg max^ 9i 
if besttheta < lasttheta then 

break the loop, 
else 

lasttheta = besttheta 
m = goodmode 
if m; = then 

AllowedSet = AllowedSet\mter/erers(/) 
end if 

mi = nextlevel(mi) 
goodmode = m 
end if 
end while 

Declare goodmode as the solution to sub problem 



This paper is submitted to NCC-2012 conference proceedings. 



4 



Draft 1 



Since we solve the problem ([T]i with a heuristic solution 
at the sub problem level of column generation, we call this 
method as Heuristic Column Generation (HCG). It should be 
noted that, this solution is sub-optimal for the sub-problem 
and hence the final solution that we obtain by using this is 
sub-optimal. 

Further, it can be predicted and also observed in sim- 
ulations presented in next section that, the final solution 
(maximum fairness) varies over a range of values, depending 
upon the initially selected random set of modes. This issue 
was not discussed much in the earlier papers using this HCG. 
Hence we in this paper, propose to choose the best value 
among solutions achieved by solving the problem for multiple 
number of initial points, analogous to the traditional tabu 
search methods used for non-convex optimization. 

In next chapter, we solve the problem for several networks 
with multiple antennas and get the final solutions. We first 
show for appropriate networks that the HCG gives a sub- 
optimal solution but is very close to the optimal solution. 
Then we see that we can solve the problem upto 30 nodes, 
almost similar to the single antennas case, as claimed in earlier 
papers. Further, we show that the MIMO does improve the 
final performance of the network and further that the gain is 
linear with the number of antennas(a) on nodes. 

V. SIMULATION RESULTS 

Consider a network of 15 nodes as shown in fig-|2] Consider 




Fig. 2. A network of 15 nodes 

the following network parameters: 
Number of Nodes = 15 
Number of Directed Links L — 
Number of antennas on each node a = 5 
Power levels used by each node = [0,4]; 
Channels are generated as L number of 5 x 5 matrices of unit 
variance, circularly Gaussian complex random numbers. 
source,destination pairs(flows) are = (7,13), (10,5), (11,8) 
The rate demands from the flows = [10, 15, 20]units 
Nodes Avg power availability = 3 units, uniform to all. 
(Note: In this scenario, we have approximately 0.5 million 
modes allowed in the network) 

By formulating the problem ([TJ and solving it directly (op- 
timally), we get the following solutions. Allowed user rates 



after solving the problem are: 

n = 5.971 units, r2 = 8.957 units 
rs = 11.942 units 
Hence the optimum fairness that can be achieved is: 
A = 0.5971 

A dominant mode which is given fairly good amount of 
time fraction is shown in fig-[3]with numbered bold arrows. 




Fig. 3. A dominant mode in the network to achieve optimal throughput 

Now we use Heuristic algorithm to solve the problem and 
see the relative performance. This is shown in fig-|4] Note that 
we use multiple initial points, we get the following solutions. 
Allowed user rates after solving the problem are: 

ri — 5.8 units, r2 ~ 8.7 units 
7-3 = 11.6 units 
Hence the optimum fairness that can be achieved is: 
A = 0.58 
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Fig. 4. Fairness Values achieved using Heuristic Algm, Compared with 
Optimal for the network of 15 nodes 

Next, we consider a network of 30 nodes and the following 
parameters. 

Number of Nodes TV = 30 
Number of Directed Links L 110 
Number of antennas on each node a = 4 
Source,destination pairs(flows) are: 
= (2,30), (5,10), (1,30), (4,6) 

The rate demands from the flows = [10, 50, 20, 60] units 
Nodes Avg power availability = 3 units, uniform to all. 
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Fig. 5. A network of 30 nodes 



Formulating the problem ([T| and solving it using Heuristic 
for multiple initial points, we see the following results fig-|6] 
Note that, this problem cannot be solved directly and hence 
the optimal performance achievable is not known. 
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Fig. 6. Heuristic algoritlimic performance for network of N=30 

Allowed user rates after solving the problem are: 
n = 2.7069 units, ra = 13.5345 units 
rg = 5.4138 units, n = 16.2414 units 

Hence the optimum fairness that can be achieved is: 
A = 0.2707 

So we demonstrated that proposed system model with JRSP 
problem is solvable upto reasonable size of the network 
and that the extended heuristic solution gives close enough 
solutions to the optimal values. We now show the performance 
gain achieved using multiple antennas in the system. 

We consider the same 15 node network and the problem 
shown in fig-|2] We vary the number of antennas on each 
node(a) from one to ten and see how the optimal performance 
(fairness) varies. Note that the channel matrices observe a 
change in dimensions as well as values, for each case. We 
then consider the same network problem but solve it using 
the Heuristic column generation. We vary the number of 
antennas (a) from one to ten and see the performance (fairness) 
variation. Both these are shown in fig-|7] We see that both 
optimal as well as Heuristic algorithmic performances improve 
linearly as we increase the number of antennas on nodes. 
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Fig. 7. MIMO gain in Optimal and Heuristic algoritlimic performance of a 
network with 15 nodes 



VI. Conclusion 

We have proposed a feasible system model for multihop 
wireless networks with multiple antennas extended from the 
MHWN with single antennas. We have shown that the problem 
is solvable upto reasonale size of the network almost similar 
to the network with single antennas. We then have shown that 
the performance gain that we achieve using multiple number 
of antennas on nodes is linear 
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